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Abstract. The observed persistence of specific periodicities detected in time 
series associated with solar surface magnetic activity over several solar cycles has 
led to numerous papers supporting the existence of preferred longitudes. Recent 
analysis of the past 120 years of sunspot number data showed that no observed 
periodicity remained coherent for durations greater than two 11-year solar cy- 
cles. Here we address the question: Could the observed periodicities of solar 
qq ■ magnetic signals on time scales of two decades be the result of a purely stochas- 

tic process? We begin to answer this by comparing phase coherence between 
observed periodic signals and signals from a model using longitudinally random 
eruptions. A surprisingly non-negligible likelihood is found, approximately 1 in 
\ 3, that observed periodicities from integrated full-disk solar parameters are a 

t^~- ■ chance occurrence for time series on the order of 20 years in duration. 
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Q ■ 1. Introduction 

The observed solar surface is spatially and temporally variable due to the dy- 
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namic nature of magnetic active regions along with differential rotation and 
meridional flows. Even with such dynamics, active regions are not observed to 
be distributed r andomly but rather magnetic regions tend to form "complexes 
of activity" (e.g. iBumba Howard! Il965l ; iBrouwer" Zwaadfl99(il h Using Kitt 



Peak Vacuum Telescope (KPV T) magnetic synoptic m aps, persistent bands of 
active "nests" were reported by Ide Toma et all ( 20001 ) to support the idea of a 



stable pattern in the tachocline for sources of emerging flux. The lifetimes of 
these activity bands were found to be as long seven solar rotations. 

For longer time scales, the temporal coherence of magnetic regions is ex- 
pected to decrease for periods greater than the lifetime of the active regions. 
However, there have been numerous reports of persi stent signals from spectra l 
analysis with the sunspot number time series (e.g. ISvalgaard Wilcox! Il975l ; 



Bogart 19821 ). U sing aut o correl ation analysis with 128 years of sunspot number 



time series data, iBogartl (19821) found a pers istent 27.5-day period signal. Ad 



ditionally, work by Neugebauer et all ( 20001 ) reported a very persistent signal 



with a well-defined period of 27.03 ± 0.002 days, using over 30 years of solar 
magnetic field and solar wind measurements. iHennev Harvevl (|2002l ) also de- 
tected this signal using photospheric magnetograms (see Figure 1). In addition, 
using phase coherence a nalysis on the past 120 years of sunspot number data, 
Hennev Harvey ( 20021 ) showed that no observed periodicity remained coher- 
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ent for durations greater than two solar cycles. Here we estimate the likelihood 
that coherent signals for durations on the order of two solar cycles are a result of 
chance occurrence. The model developed and used for this study can operate in 
several modes ranging from purely random longitudinal location for the infusion 
of new flux to including a specific longitude and rotation rate for new flux at 
any selected latitude. Since the integrated full-disk signal is of interest here, the 
solar model was used in its random l ongitude mode and does not include de- 
tailed solar surface flux transport (e.g. Schriiver fc DeRosa 120031 ) . The following 
section outlines the magnetic activity model utilized for this preliminary study. 



2. Magnetic Activity Model 

The modeling of large-scale solar magnetic activity is achieved in several steps. 
The initial step is to estimate the values of the mean eruption strength and decay- 
time for observed flux within a given latitude band. These values were deter- 
mined using the same unsigned magnetic synoptic ma ps created from KPVT ful l- 
disk photospheric measurements that were us ed by Henney &: Harvey ( 2002i ). 



The KPVT synoptic map data is described bv lHarvev et al. ( 1980l ). The long- 



term trends of solar magnetic activity are determined for 18 latitude bins of 
the KPVT synoptic magnetic maps. The maps were spatially sampled in sine- 
latitude bins and temporally sampled to create latitudinal time series. The 
temporal gaps are filled with a high-order autoregressive model, and each time 
series is temporally filtered using a non-recursive digital filter with a low-band 
pass between 10 and 60 days. Each time series is then median filtered with a 
window size of 83-days, then finally smoothed by convolving a Bartlet function 
with a window size of 27-days. 

The filtered-flux time series is sampled at a cadence of 10 days, then is 
interpolated to a final cadence of 1 day. Allowing for daily eruptions, the model 
eruption strength and decay-time for each latitude band are then estimated by 
minimizing the standard deviation between the filtered- and modeled-flux time 
series (see Figure 2). The eruption strengths shown in Figure 2 vary from ap- 
proximately 1 x 10 22 to 3.5 x 10 22 Mx. These values are in agreement with 
previous estimates of 1 x 10 22 Mx and 2 x 10 22 Mx for observed ephemeral (e.g. 
lHarvev Sz Martenlll979l ) and active (e.g. ISheelevlll966l ) region strengths respec- 



tively. The estimated decay periods deviate around 100 days, which correspond 
to approximately 1 x 10 20 to 3 x 10 20 Mx/day. These de cay rates are c onsistent 
with previous estimates of 2 x 10 20 Mx/day reported by iHarvevI (1993). 



The modeled flux emergence is governed by the input filtered-flux time se- 
ries. After an initial period of approximately one hundred days of flux emergence 
to insure that equilibrium has been reached with the filtered-flux, the modeled 
time series is created for the same interval and cadence as the observed data set 
used to create Figure 1. Each latitude bin is assumed to rotate in agreement 
with the rotational profile observed by iHoward &: Harvevl ( 19701 ). New flux is 



injected into a given latitude bin if the model flux drops below the filtered-flux 
time series for that latitude. The eruptions occur randomly in longitude over a 
range between and 2tt. The eruption strengths are random with a distribution 
and mean value estimated from the filtered-flux time series. After a period of 
five times the decay-time the contribution from any given flux emergence is set 
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Figure 1 . Comparison of power and phase coherence spectra between Kitt 
Peak Vacuum Telescope (KPVT) unsigned magnetic flux integrated full-disk 
time series (top t wo), and the internation al sunspot number (bottom two) 
time series (from iHennev fc Harvevl I2002T ) . Each time series is sampled at 
a cadence of 1 day and for the same 8839-day interval (1977.2 to 2001.4, in 
units of fractional years). The 3-er noise level for the phase coherence spectra 
is delineated by the horizontal dashed lines. The vertical dotted line illustrates 
the 27.03-day period position. 
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Figure 2. Emergence strength (top), decay time (middle), and number of 
eruptions (bottom) per latitude bin for the magnetic activity model used to 
create the simulated time series. These values were estimated by minimizing 
the standard deviation between the filtered- and modeled-flux time series. 



to zero. For the 8839-day interval simulated, there were approximately 4.84 x 10 4 
eruptions distributed over the 18 latitude bins for each model realization. 



3. Modeled Time Series Analysis 

The solar magnetic activity model was run for 500 realizations, each case creat- 
ing a simulated integrated full-disk magnetic flux time series with a total of 8839 
days at a cadence of 1 day. Each time series was tested for signals that show sim- 
ilar levels of coherence to those illustrated in Figure 1, where a fixed frequency 
harmonic signal is considered coherent if it has negli gible temporal variance i n 
phase for the interval of time investigated. Following iHennev Sz Harvevl (2002), 



we estimate the phase and amplitude of a signal from a non-linear least-squares 
fit of a fixed frequency cosine function for overlapping segments of a time se- 
ries. The phase can also be obtained with Fourier analysis, with the caveat 
that the phase error est imate is not readily available (for more discussion, see 
Hennev Harvey 20021 ). Examples of the power and phase coherence spec- 



tra for the simulated unsigned magnetic flux integrated full-disk time series are 
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Figure 3. Comparison of power (top and second from the bottom) and 
phase coherence (bottom and second from the top) spectra between two sam- 
ple simulated unsigned magnetic flux integrated full-disk time series. Time 
series with phase coherence values of approximately 0.02 (top two) and 0.04 
(bottom two) respectively are exhibited. The 3-ct noise level for the phase co- 
herence spectra is delineated by the horizontal dashed lines. Each simulated 
time series is for the same interval and cadence as in Figure 1. 
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Figure 4. Each of the 500 realizations of the modeled solar magnetic activity 
time series, the signal having the maximum phase coherence is selected. For 
these signals the distribution of signal-to-background (top), phase coherency 
(middle), and period (bottom) of the maximum phase coherence spectrum 
signal are shown. The vertical dotted line illustrates a signal-to-background 
ratio of 10 (top) and a phase coherence of 0.012 (middle). See text for dis- 
cussion. 



shown in Figure 3. The top two and bottom two figures, from simulated time 
series with phase coherence values of approximately 0.02 and 0.04 respectively, 
are similar to those shown in Figure 1. 
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4. Results and Discussion 

As a measure of the similarity between the simulated time series and those 
observed, the signal having the maximum phase coherence is selected for each 
of the 500 realizations. The distribution of the measured signal-to-background 
(S/B), the phase coherency, and period of the maximum phase coherence spec- 
trum signal from eac h realization of the mode led solar magnetic activity are 
shown in Figure 4. In lHennev Harvevl (2002). it was shown that a phase 



co- 



herence value greater than 0.012 or with a S/B greater than 10 for the length 
of time investigated here corresponds to a signal that is phase coherent for ap- 
proximately over 50% of the interval investigated. Illustrated by the vertical 
dotted lines in Figure 4, roughly 33% of the maximum phase coherent signals 
have a S/B greater than 10. In addition, nearly 34% of the signals have a phase 
coherence value greater than 0.012. By comparison, the 27.03-day signal in the 
sunspot number time series in Figure 1 has phase coherence and S/B values of 
approximately 0.016 and 7 respectively. 

If the solar magnetic model also included act i yity " nesting" on time scales 
up to 7 months, as reported by Ide Toma et al.l ( 20001 ). we would expect co- 



herence values approaching those of the observed magnetic time series. This 
is supported with an extreme test using the current model in a mode that al- 
lows for coherent eruptions at a fixed period of rotation for a chosen latitude, 
while still matching the observed mean flux for each latitude bin as discussed in 
Section 2. This resulted in notably higher than observed peaks at the selected 
rotational period in both the power and coherence spectra. Observationally, the 
27.03-day signal shown in Figure 1 is found to be phase coherent for durations 
of three years for terms scattere d throughout the 24-year time series (see Fig- 
ure 2 in Hennev Harvey 2002). This suggests that activity nesting may be 



sustained to form long lived magnetic complexes (LLMC) lasting for periods of 
a few years. This tentatively suggests that if the LLMC are the explanation for 
the observed coherence, they are scattered throughout the entire time series. If 
this is the case, the LLMC must rotate with the same period to contribute to 
the coherence spectrum. Tentative arguments favoring this suggestion can in 
fact be advanced: LLMC could resist dispersion if they have not broken their 
connection to the lower solar convection zone. Conceivably then, the generation 
of such fields can be favored at a given latitude. Work aimed at clarifying these 
issues is in progress. 



5. Conclusion 

We have shown that random eruptions of magnetic flux, modulated with solar 
differential rotation, have a surprisingly high probability of resulting in a coher- 
ent signal comparable to observation. The solar magnetic model utilized here 
with longitudinally random eruptions is treated as a baseline case since addi- 
tional processes must be included to explain observed nesting of magnetic activ- 



i ty on time scales of months to a few years. For example, in iHennev Harvey 



(2002) it was found that the origin of the 27.03-day signal is most likely the result 



of long-lived, on the order of a few years, complexes of magnetic activity in the 
solar northern hemisphere. Whether long-lived complexes of magnetic activity 
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merely by themselves are sufficient to explain the observed periodicities is under 
study . The inclusion of activity complexes or nesting, e.g. ISchriiver Sz DeRosal 
(2003), with variable lifetimes and latitudes is under study and will be presented 
in future work. With the inclusion of nesting processes in the model, it is ex- 
pected that the probability of resulting coherent signals would increase. This 
strongly suggests that spectral analysis results with integrated full-disk param- 
eter signals, e.g. the sunspot number time series, should be interpreted with 
care. It should be stressed that these results do not lessen the possibility of 
very long-lived (between 30 to 300 solar rotations) magnetically active longi- 
tudes. However, these findings reinforce the need for spatial source tracking and 
phase coherence analysis to better understand the source of these signals and 
substantiate claims of very long-lived active longitudes resulting from spatially 
integrated time analysis of solar parameters. 
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